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In this talk, we report on a recent next-to-leading order QCD calculation of the production of 
a W boson in association with three jets at hadron colliders. The computation is performed 
by combining two programs, BlackHat for the computation of the virtual one-loop matrix 
elements and Sherpa for the real emission part. The addition of NLO corrections greatly 
reduces the factorization and renormalization scale dependence of the theory prediction for 
this process. This result demonstrates the applicability of unitarity-based methods for hadron 
collider physics. 

1 Introduction 

The production of a vector boson in association with jets is an important process at the LHC. 
Apart from its interest as a test of QCD, it contributes significantly to the background of many 
Standard Model processes {tt production, single top production, and Higgs decay to two vector 
bosons) and new physics processes. Successful measurements of these processes require a reliable 
theoretical description of the vector boson -|- jets processes. 

Leading-order QCD predictions for processes with jets suffer from a large dependence on the 
renormalization and factorization scales. This problem can be tamed by adding next-to-leading 
order (NLO) corrections. Such corrections are composed of two parts. The real corrections to 
an n-parton process arise when an additional parton is emitted, in an (n -|- l)-parton process. 
One-loop n-parton amplitudes generate the virtual part of the NLO correction. 

BlackHat ^ is a numerical implementation in C-|--|- of so-called on-shell methods for comput- 
ing one-loop amplitudes. The starting point for a one-loop amplitude A with massless propaga- 
tors is its general decomposition in terms of scalar integrals, 

A = J2c^Il + J2clIl + Y.clP, + R , (1) 

i i i 

where Z^, /I, II are scalar bubble, triangle and box integrals. The 'rational term' i? is a rational 
function of spinor products and does not contain any logarithms. The objective of on-shell 
methods is to determine the coefficients of the integrals and R without using Feynman 
diagrams. We refer the reader to the literature^ for more details on these methods. 

In our numerical implementation, the coefficients of the integrals are determined using the 
analytic approach of Forde^, which is related to other recent methods^. The rational term is 




Figure 1: NLO Diflterential cross section do-{W — > ei/+ > n — jets) /dEJ^^^^^^^ for n = 1,2 compared with the 
measured cross section. The upper panels contain the LO and NLO parton level predictions and the CDF data 
points with their statistical and total uncertainties represented by the inner and outer error bars respectively. 
The distributions normalized by the NLO prediction are shown in the lower panels. The scale uncertainty of the 
different predictions is represented by the shaded grey (NLO) and orange (LO) bands. The dotted black line 

represents our leading color approximation. 
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Table 1: Comparison of the total cross sections in pb for W + n jets with _Ey* > 25 GeV from CDF to 
NLO QCD. For 1 and 2 jets the cross sections with and without LC approximation are displayed to show the 
quality of the approximation. For the three jets result, only the LC NLO result is currently available, but we 
expect a similarly small deviation for the full NLO result. The experimental statistical, systematic and luminosity 

uncertainties have been combined for the CDF results. 



computed using on-shell recursion relations for one-loop amplitudes^. Numerical stability of the 
implementation is achieved by using high-precision libraries ^ when (and only when) necessary. 
This stability has been demonstrated elsewhere'^. 

Sherpa*^ is a C++ Monte Carlo event generator. It can compute the real part of the NLO 
corrections in an automated way ^ using Catani and Seymour's dipole subtraction method 
In addition, the subtraction term integrated over the unresolved phase space is provided. The 
results presented below have been integrated over the relevant phase space using Sherpa. 

For the + 3 jets virtual cross section we used a leading-color (LC) approximation for the 
finite part of the virtual amplitude. This approximation amounts to neglecting the terms in the 
ratio of the virtual terms to the tree cross section that are suppressed by factors of ^/N^ (color 
suppressed) or Nf/Nc (virtual quark loop). This approximation has been shown to be very good 
in the following section for W + 1,2 jets , so we expect it to be valid for + 3 jets. A related, 
but different, approximation that includes only a subset of partonic subprocesses has been used 
in another computation^^ of 14^ + 3 jets. The benefit of our approximation is that the number of 
(color ordered) primitive amplitudes to evaluate is significantly reduced. We checked agreement 
between the primitive amplitudes we used here and those found in a different calculation . 



2 jets at the Tevatron 

We compare the NLO prediction for W + 1,2,3 jets with data from the CDF experiment ^'^ at 
the Tevatron. For the analysis we have used the same cuts as in the CDF analysis with the 
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Figure 2: The left panel shows the NLO Differential cross section da{W -+ ei^+ > 3 — jets)/dE^'^~''^^ compared 
with the measured cross section. Its upper part shows the LO and NLO parton level predictions (within the 
LC approximation for the virtual part) and the CDF data points with their statistical and total uncertainties 
represented by the inner and outer error bars respectively. The normalized distributions are shown in the lower 
part of the left panel. The scale uncertainties are represented by the shaded grey (NLO) and orange (LO) bands. 
The right panel shows the scale dependence of the total cross section on the renormalization and factorization 
scale jj,, taken equal and varied between 1/4 and 4 times the mass of the W boson. 



SISCone jet algorithm instead of the JETCLU cone algorithm used by CDF, as the latter 
is not infrared safe. We set an event-by-event renormalization and factorization scale according 
tofi = ^mlr + P^{W). 

In Figure 1 we present the transverse energy distribution of the n-th jet for W + 1,2 jets 
jets production. The NLO result agrees with the previously available results from MCFM ^'^ . 
The lower part of these plots shows the reduced scale dependence of the NLO prediction. The 
dotted black line demonstrates the validity of our LC approximation across the whole Ej- range. 
The plot in Figure 2 shows a good agreement between the NLO prediction (within our LC 
approximation) for the Et distribution of the third jet and the experimental data. The second 
plot of Figure 2 displays the improvement of the scale dependence of the cross section when 
NLO corrections are added. More details on the setup of our analysis can be found elsewhere 



3 W+3 jets at the LHC 



We repeated the same analysis for the LHC with a center-of-mass energy of 14 TeV. For this 
analysis, we chose cuts suggested by the ATLAS and CMS technical design reports: Ej^ > 20 
GeV, \rj^\ < 2.5, fir > 30 GeV, M]f > 20 GeV, and E^^^ > 30 GeV. Here the E'^s are 
transverse energies, is the missing transverse energy, the transverse mass of the eu 

pair, and rj the pseudorapidity. The -ET-ordered jets are required to have a rapidity in the range 
\r]\ < 3. We used SISCone with R = 0.4. Figure 3 shows the distributions in the scalar 

transverse energy sum Ht and in the three-jet mass Mjjj = {kji + kj2 + fejs)^ for W~ + 3 jets 
at the LHC. 
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Figure 3: The NLO prediction for the scalar transverse energy sum Ht (left panel) and the three-jet mass Mjjj 
(right panel) compared to the LO prediction. The LC approximation has been used for the virtual contribution 
to the NLO result. The scale uncertainties are represented by the bands in the lower parts of the plots. 
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